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Abstract 



Using the South Pole Acoustic Test Setup (SPATS) and a retrievable transmitter deployed in holes drilled for the 
IceCube experiment, we have measured the attenuation of acoustic signals by South Pole ice at depths between 190 
m and 500 m. Three data sets, using different acoustic sources, have been analyzed and give consistent results. The 
method with the smallest systematic uncertainties yields an amplitude attenuation coefficient a — 3.20 ± 0.57 km -1 
between 10 and 30 kHz, considerably larger than previous theoretical estimates. Expressed as an attenuation length, 
the analyses give a consistent result for A = 1/q of ~300 m with 20 % uncertainty. No significant depth or frequency 
dependence has been found. 

Key words: neutrino astronomy, acoustics, South Pole, acoustic attenuation, ice 
PACS: 47.35.Rs, 62.65. +k, 92.40.Vq, 93.30.Ca, 95.55.Vj, 92.40.vv, 91.60.Qr, 96.25.hf 



1. Introduction 

Experiments to study ultra- high-energy neutrinos 
have been the subject of increasing interest during 
recent years [1-3] . As observed by HiRcs [4] and the 
Pierre Auger Observatory [5], the charged cosmic 
ray flux decreases steeply above 10 19,5 eV. This is 
most likely due to the interactions of charged cosmic 
particles with the cosmic microwave background ra- 
diation, known as the GZK effect [6,7]. The detec- 
tion of neutrinos from this interaction would con- 
firm this explanation. Spectral, temporal, and direc- 
tional distributions of such neutrinos, enabled by de- 
tecting a significant number of them, would address 
important questions of cosmology, astrophysics and 
particle physics [8] . 

Estimates of the small flux of GZK neutrinos vary 
by an order of magnitude. The results from Engel, 
Seckel, and Stanev [9] are often used MS (1 standard 
for the discussion of possible detector scenarios. In 
all cases, detector effective volumes of 100 km 3 or 
larger are required. 

Several past and current experimental projects 
seek to detect GZK neutrinos. The best limits cur- 
rently come from searches for radio signals from 
neutrino interactions. Presently the most stringent 
flux limit is from the ANITA project [10]. However, 
all experiments searching for weak particle fluxes 
must contend with challenging systematic effects 
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and background separation. This may be overcome 
in the future by using a hybrid detector. The radio 
technology could be complemented by adding acous- 
tic detectors, searching for the sound produced by 
the interactions of neutrinos above 10 18 eV [11,12]. 
Simulations of a radio-acoustic hybrid detector en- 
compassing the optical IceCube neutrino observa- 
tory at the South Pole [13] , assuming an attenuation 
length on the order of a few kilometers as theorized 
in [14], gave a predicted detection rate of 20 neu- 
trino events per year, half of them detected by both 
radio and acoustic sensors [15]. 

Ice seems to be a favorable medium for the ap- 
plication of the optical, radio and acoustic detec- 
tion techniques. Ice properties have been measured 
for optical signals by the AMANDA experiment [16] 
and for radio waves by the RICE experimental pro- 
gram [17,18]. As far as the propagation of acoustic 
signals is concerned, theoretical estimates [14] have 
indicated low absorption and scattering correspond- 
ing to an attenuation length greater than 1 km. 

To test the theoretical estimates, the South Pole 
Acoustic Test Setup (SPATS) was deployed in Jan- 
uary and December 2007 in the shallowest 500 m of 
IceCube holes at the South Pole. SPATS is used to 
study the following four properties: 

- the speed of sound as a function of depth in the 
ice, important for the expected neutrino signal 
strength, event localization and reconstruction; 

- the acoustic noise level at the South Pole, deter- 
mining the energy threshold of a future neutrino 
detector; 

- the rate and nature of transient acoustic signals 
which could mimic neutrino interactions, present- 
ing a serious background source; 

- the attenuation of acoustic signals by the ice, 
which determines the necessary density of acous- 
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tic sensors in the ice for a reasonable detection 

efficiency for neutrino interactions. 
The sound speed has been measured in [19]. Noise 
and transients are currently under study. Our mea- 
surement of attenuation by the ice is reported here. 

2. Experimental setup 

2.1. The South Pole Acoustic Test Setup 

Figure 1 shows the geometry of the acoustic in- 
strumentation used for the measurement. In this 
section we provide a description of all the hardware 
used, while in the next we explain which sound 
source and which subset of sensors we used for each 
measurement. The permanently deployed in-ice 
SPATS hardware consists of four vertical instru- 
mented cables (strings), which were deployed in the 
shallowest 500 m of four IceCube holes. The esti- 
mated surveying error on the position of each hole 
center is ±0.1 m in each of the x and y coordinates. 
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Fig. 1. SPATS geometry and IceCube footprint in 2009. 
The permanently deployed SPATS hardware is indicated by 
open circles with the string ID (ABCD) and corresponding 
IceCube hole number. The full circles indicate the location 
of the holes, indexed by the corresponding IceCube hole 
number, where the pinger was deployed in the 2008-2009 
season. 

The maximum string-to-string horizontal dis- 
tance is 543 m between String C and String D. 
Strings A, B and C have acoustic stages (described 
below) installed at 80, 100, 140, 190, 250, 320 and 
400 m depth. String D has instrumentation at 140, 



190, 250, 320, 400, 430 and 500 m depth. A diagram 
of the deployed instrumentation is presented in [19]. 

An acoustic stage consists of a transmitter mod- 
ule and a sensor module. The transmitter module 
consists of a steel pressure vessel of ~10.2 cm di- 
ameter that houses a high- voltage (HV) pulse gen- 
erator board and a temperature or pressure sensor. 
The active element, a ring-shaped piezo-ceramic el- 
ement cast in cpoxy for electrical insulation, is po- 
sitioned ~13 cm below the steel housing. The HV 
pulse stimulates the piezo-ceramic element, result- 
ing in an acoustic pulse. The amplitude of the acous- 
tic pulse can be modified externally by adjusting two 
additional input signals: the steering voltage and the 
trigger pulse length. 

The SPATS sensor module has three channels, 
separated by 120° in azimuth to ensure good an- 
gular coverage. A channel consists of a cylindrical 
piezo-ceramic element that is directly soldered to a 
three-stage amplifier and pressed against the steel 
housing. The three channels of each sensor are inde- 
pendent and therefore can be treated as three stand- 
alone receivers. From previous laboratory tests, it is 
known that the sensitivity of each sensor depends 
on both the polar and the azimuthal angles [20]. 
Most of the sensors were calibrated in water before 
installation, but there is evidence that the sensitiv- 
ities have changed since deployment. The effect on 
the sensor sensitivity of simultaneous low tempera- 
ture, high pressure and ice coupling (compared to 
water coupling where calibration was performed in 
the laboratory) is unknown, as is the precise loca- 
tion of the ~10.2 cm diameter sensor vessels in the 
75 cm wide IceCube hole. During re-frcczing of the 
holes, ice with bubbles and perhaps cracks is formed 
around the sensors, possibly leading to a strong but 
unknown angular dependence. A method to cali- 
brate the sensors' beam pattern in situ has not yet 
been developed. 

The Hydrophone for Acoustic Detection at the 
South Pole (HADES) was designed and constructed 
to offer an alternative sensor with a different dy- 
namic range. A ring-shaped piezo-ceramic element 
is connected to a two-stage differential amplifier 
placed inside the ring. The assembly is located out- 
side of the housing of the SPATS sensor module and 
is coated with polyurethane plastic resin in order 
to protect the electronics from water and ice. Two 
HADES modules, each with one sensor channel, are 
installed in place of SPATS sensors on String D, at 
190 m and 430 m depth. 

We refer to each sensor channel by the string iden- 
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tifier letter (ABCD), the number of the stage (1-7, 
from shallowest to deepest) and the number of the 
channel (0-2). For example, AS6-0 indicates channel 
of sensor module number 6 of String A. 

Each of the acoustic stages is connected to an 
Acoustic Junction Box (AJB) at the surface. This 
AJB is a strong aluminum box buried under ^2 m 
of snow and contains a rugged embedded computer 
(String PC) with the electronic components nec- 
essary to digitize locally all signals. Each AJB is 
connected to a computer (the Master PC) , located 
in the IceCube Laboratory, which stores the data 
until it is transferred by satellite to the Northern 
hemisphere or copied locally to tape. A GPS-based 
LRIG-B time code signal provides absolute time 
stamping. For a detailed technical description of the 
SPATS permanent hardware and its development, 
see [19,20]. 

2.2. Retrievable pinger 

In addition to the equipment deployed in the ice, 
a retrievable transmitter (pinger) was designed as 
a unique acoustic source to be operated in multiple 
water-filled IceCube holes, prior to IceCube string 
deployment. A previous version of the pinger was 
used in the 2007-2008 season to measure the sound 
speed profile. Here we describe the pinger as mod- 
ified for the attenuation measurement in the 2008- 
2009 season, when it was deployed in IceCube holes 
5, 19, 28 and 37 (Figure 1). 

The pinger is an autonomous transportable de- 
vice consisting of an acoustic stage which is lowered 
into the water, and an auxiliary box, sitting on the 
surface, which provides power and a trigger signal. 
The box was connected to the acoustic stage during 
operation through a steel-armored four-conductor 
cable, about 2700 m long, spooled on a winch, used 
to lower and raise the pinger. 

The acoustic stage consists of a custom designed 
high- voltage pulser circuit and of a spherical piezo- 
electric ceramic emitter 4 from the International 
Transducer Company (ITC). The Transmitting 
Voltage Response (TVR) specification provided by 
ITC has dominant components in the frequency 
range between 10 and 30 kHz, with a peak of 
149 dB re (/zPa/V @ 1 m) at the resonance peak. 
The resonance frequency has been measured to be 
f res = 17.7 kHz. 
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The box on the surface, (the Acoustic Pinger Box, 
APB) contains a 24 V scaled lead acid rechargeable 
battery pack (Hawker-Cyclon) rated to low temper- 
atures. A GPS clock (Garmin model GPS 18 LVC) 
generates a 1 pulse-per-second (PPS) signal. A fre- 
quency multiplier circuit, consisting of a Complex 
Programmable Logic Device (CPLD) and an oscilla- 
tor, produces a continuous train of electric pulses at 
a rate which is adjustable by a switch. The rate of the 
pulses was 10 Hz (in holes 28, 19 5) or 8 Hz (in hole 
37) during operation at the South Pole. This means 
that a sequence of pulses equally spaced in time (by 
100 or 125 ms) was produced by the CPLD, with ev- 
ery 10th (or 8th) pulse synchronized with the PPS 
rising edge. Each electric pulse is delivered through 
the APB to the acoustic stage. Here it triggers a 
timer in monostable configuration which charges an 
LC circuit for a defined time, at the end of which 
the energy accumulated is transferred to the piezo- 
electric ceramic, resulting in acoustic emission. The 
electric pulse exciting the piezoelectric ceramic is 
unipolar, about 300 V high, and has a full-width- 
half- maximum width of 60 /is. The pressure signal 
transmitted to the ice is the convolution of the stim- 
ulating electrical pulse and the TVR. The calcula- 
tion of the spectrum, displayed in Figure 2, shows 
that most of the power is emitted in two lobes, one 
around 10 kHz and one around 20 kHz, with a mini- 
mum at the resonance frequency of the piezoelectric 
ceramic. 




Frequency (kHz) 



Fig. 2. Spectrum of the acoustic pinger pulse, calculated as 
the convolution of the electrical high-voltage pulse and the 
Transmitting Voltage Response of the ITC-1001. 

The pinger was used in 2007-2008 to measure the 
sound speed vs. depth in South Pole ice for both lon- 
gitudinal or pressure waves (P waves) and transver- 
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sal or shear waves (S waves) [19]. The production 
of shear waves (from mode conversion at the water- 
ice boundary) in that data set was favored by the 
fact that the pingcr's lateral position in the hole was 
off-center and varying. However, this made atten- 
uation measurement difficult due to the change in 
shape of the waveform recorded by one sensor chan- 
nel for different pinger holes. For the 2008-2009 sea- 
son, the acoustic stage was equipped with mechani- 
cal centralizers, suitable to keep the acoustic emitter 
close to the central axis of the hole (see description 
in [21]). This prevented the stage from swinging and 
stabilized the acoustic pulse transmitted in the ice. 
The centralization of the pinger (and the generally 
longer distances between pinger and sensors in 2008- 
2009 than 2007-2008) minimized the appearance of 
shear waves in the 2008-2009 data set; consequently 
it has been possible to measure the attenuation of P 
waves without the complication of S waves. 

3. Data processing and analysis techniques 

3.1. Data sets 

Three different types of data were acquired and 
analyzed to measure the acoustic attenuation: 

(i) Pinger data: multiple acoustic pulses were 
emitted from several water-filled IceCube 
holes and recorded by sensors of the SPATS 
array. For each pinger depth, 180 (144 in hole 
37) single acoustic pulses were recorded by 
each sensor channel. 

(ii) Inter-string data: transmitters of the SPATS 
array emitted acoustic pulses that were 
recorded by the SPATS sensors. For each 
transmitter 500 pulses were recorded by each 
single sensor channel. 

(iii) Transient data: sound produced in re-freezing 
IceCube holes at depths of about 250 m was 
recorded by the sensors of the SPATS array. 

An overview of hardware components involved in the 
analysis, together with the corresponding range of 
distances (baselines) and number of pulses recorded 
for each channel, is shown in Table 1. 

3.2. Data processing 

For each of the analyses described below, N p con- 
secutive acoustic pulses from one source (transmit- 
ter or pinger) are recorded by one sensor channel 



Table 1 

Overview of data used for the attenuation measurement. 
The last column gives the number of pulses in each sensor 
recording. 



Data type Receiver 


Source 


Distance 


Number of pulses 


Pinger SPATS 
sensors 


pinger 


125-1023 m 


144-180 


Inter-string SPATS 
sensors 


SPATS 
transmitters 


125-686 m 


500 


Transients SPATS 
sensors 


transient 
events 


243-750 m 


1 



and are averaged to improve the signal-to-noise ra- 
tio. Each string uses a single clock to drive both its 
analog-to-digital converters (ADCs) and its digital- 
to-analog converters (DACs). The ADCs are used 
to record the sensor waveforms with a sampling fre- 
quency of 200 kHz. The DACs are used to pulse 
SPATS transmitters. The clocks drift at a rate that 
is typically several parts per million, or tens of /xs 
over the duration of a single sensor recording. This 
cumulative amount of drift is on the order of one 
signal oscillation period and therefore can cause se- 
vere decoherence in pulse averaging if the nominal 
rather than true sampling frequency is used. On the 
sensor side we correct the clock drift effect by us- 
ing the IRIG-B GPS signal (which is recorded syn- 
chronously with each sensor channel) to determine 
the actual sampling frequency at the time of the 
recording. We then use this actual sampling fre- 
quency to average the recorded pulses. 

The pinger pulse emission is driven by the GPS re- 
ceiver and by the frequency multiplier. The GPS re- 
ceiver provides the PPS pulse with a delay of about 
1 /is relative to the GPS signal which drives the 
String PC ADCs. The frequency multiplier synchro- 
nizes the train of electric pulses with the PPS pulse, 
but introduces some jitter in time which we mea- 
sured in the laboratory to be ± 5 /is over the recorded 
duration for a single channel. Time jitter smears out 
the amplitudes contributing to each single point in 
the averaged sensor waveform, the effect of which is 
automatically included in the statistical uncertainty 
of the average pulse amplitude. 

Inter-string data (recorded with the frozen-in 
transmitters rather than the retrievable pinger) 
have the additional complication that the transmit- 
ters are not synchronized to GPS. These transmit- 
ters pulse at a rate which can drift relative to abso- 
lute time. Therefore in this data, in addition to the 
drift of the recording sensor, we must correct for the 
clock drift on the transmitting string, which causes 
the actual transmitter repetition rate to be different 
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from the nominal one. For each run we measured 
the actual repetition rate using IRIG-B GPS signals 
in the same way as described above for the sensors. 
The true repetition rate varied from string to string 
but was constant on each string at the few-percent 
level over the course of the two-day inter-string data 
taking campaign. Therefore for each string we used 
the mean value of the transmitter repetition rate in 
averaging the recorded pulses. The effect of residual 
clock drift due to using the average drift rate at each 
transmitting string rather than using the instanta- 
neous drift rate is < 0.3 sample over the duration of 
a 20 s sensor recording, much less than a single sig- 
nal oscillation period and therefore not sufficient to 
cause any residual decoherence in pulse averaging. 

After correcting the sample times for clock drift, 
we wrap the waveform samples in time modulo the 
pulse repetition period, in order to overlay all N p 
pulses recorded in a single sensor recording. We then 
time-order the samples and bin them in groups of 
Nbm consecutive samples. The number of samples 
per bin N^n was chosen to be one-half the number 
of pulses for the inter-string data, and equal to the 
number of pulses for the pinger data. The choice 
of bin size was made independently for the inter- 
string and pinger data because they have different 
frequency content and signal-to-noisc ratio. For each 
bin j, we compute the mean voltage Vj and the stan- 
dard error of the mean try^. , which we use as an es- 
timate of the uncertainty on the sample amplitude 
because the noise is observed to be stable and Gaus- 
sian. The resulting average pulse is one repetition 
period in duration. For a detailed description of the 
algorithm used sec [22,21]. 

3.3. Analysis techniques 

The basic quantity measured by the sensors is a 
pressure wave exciting their piezo-ceramics for a du- 
ration which depends on the sensor construction, 
the signal strength and the distance between emitter 
and receiver (see Figures 3(a) and 3(b)). The signal 
amplitude A recorded by a sensor is proportional to 
the input acoustic pressure. For a point source with 
spherical emission, A depends on the distance to the 
source and the amount of attenuation by the ice: 

A = ^er ad = ^e- d '\ (1) 
a a 

where A$ is constant for a given transmitter and 
sensor channel and depends on the amplitude of the 
sound emitted at the source and the sensitivity of 



the receiver; d is the distance to the source; a is the 
acoustic attenuation coefficient and the attenuation 
length A is its inverse. We assume here that a is 
independent of frequency, position, and direction. 
We multiply both sides of Equation 1 by the known 
distance d and take the natural logarithm to define 
a new variable y: 

y = HM (2) 

In this way we can turn the previous nonlinear equa- 
tion into a linear one: 

y — In Aq — ad — —ad + 6, (3) 

where b is a free normalization parameter related to 
the sensitivity of the particular sensor piezoceramic. 
Linear regression can be applied to determine the 
best fit and uncertainty on each of the two parame- 
ters a and b. 

The measured signal has a complex waveform that 
depends on the particular choice of transmitter and 
sensor (see Figures 3(a) and 3(b)). It is difficult to 
use one particular peak amplitude of this waveform 
directly to characterize an arriving signal. A more 
stable quantity is the energy of the recorded pulse, 

n 
i=l 

where Vi is the amplitude of the sample i (after av- 
eraging, as explained in section 3.1) and the number 
of samples n is waveform-dependent. Because every 
sample i of the mean waveform has a statistical un- 
certainty oy i , we can calculate a statistical uncer- 
tainty <je for the energy. We note that this quantity 
is an "energy" in the signal-processing sense and is 
directly proportional to the acoustic energy in the 
pressure pulse. 

From the energy E, once the noise has been sub- 
tracted as explained below, we can also calculate the 
effective amplitude of the pulse: 

A eff = VE (5) 

and thereby extract the amplitude attenuation co- 
efficient from Equation 3 (using A = A e ff) rather 
than the energy attenuation coefficient. The statis- 
tical uncertainty of the effective amplitude is de- 
termined with standard error propagation from ue- 
This method has been used for most of the studies 
described below. 

An alternative approach is to apply Eqs. 1 - 3 to a 
calculation of the waveform energy in the frequency 
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domain. In this case, the effective amplitude A e ff is 
given by 




with the coefficients of the noise-subtracted Fourier 
spectrum S m . Results of both methods will be given 
below. 

In both methods, we do not deal directly with the 
response function of the sensors, which is implicitly 
included in Ag. The frequency-dependent sensitiv- 
ity of our sensors exhibits several resonances due 
to both the piezoelectric ceramics and the housing, 
with their corresponding couplings. As discussed in 
section 4.2, in the pinger measurements we assume 
that the sensor response function is constant over 
time for all waveforms recorded with the same sen- 
sor channel (therefore Equations 1-3 can be applied 
exactly since Aq is really independent of d). In the 
other analyses (inter-string and transient data) the 
variation of Ao is included as a source of systematic 
uncertainty. 

4. Pinger measurement 
4.1. Pinger data acquisition 

As described in Section 2.2, the pinger was oper- 
ated in four IceCube holes (28, 19, 5, 37), shown in 
Figure 1, in the 2008-2009 season. The pulse repeti- 
tion rate was 10 Hz in all holes except for hole 37, for 
which the repetition rate was 8 Hz. The distances 
from the SPATS array ranged from 156 m (String 
C to Hole 37) to 1023 m (String D to Hole 5). The 
range in azimuth angles spanned by the pinger holes 
as seen from an acoustic sensor was small: 13 degrees 
for String D, 7.2 degrees for String A, 6.6 degrees 
for String B, and 8.2 degrees for String C. 

The pinger was lowered from the surface to a max- 
imum depth of 500 m and then raised back to the 
surface. Both on the way down and on the way up, 
the pinger was stopped for ^5 min at nominal depths 
of 190, 250, 320, 400, 500 m, which arc instrumented 
with SPATS sensors. The depth was monitored us- 
ing the cable payout, initially calibrated by count- 
ing the number of turns of the winch. In addition, 
multiple calibrated pressure sensors attached to the 
acoustic stage recorded the hydrostatic pressure as 
a function of time. After deployment we averaged 
the pressure and the payout data and verified that 



the stop depths were within ± 5 m of the nominal 
values. 

Only those recordings of sensors at the same depth 
of each pinger stop depth were analyzed for atten- 
uation, in order to perform the measurement along 
horizontal paths (see more details in the section 4.2) ; 
all three channels of the same sensor were recorded 
simultaneously (differently from [19]). 

Due to the time necessary to transfer the data 
from the String PC to the Master PC through a DSL 
connection over a surface cable, the optimal dura- 
tion of the recording at 200 kHz on the 3 channels 
of one sensor was found to be 18 s and this duration 
was used for the pinger runs. 

As can be seen in Figure 3(a), two waveforms 
recorded while the pinger was stopped during low- 
ering and raising, respectively, within the same 
hole look very similar and overlap very well. This is 
thanks to the mechanical stabilization of the pinger 
in the hole, in contrast to the 2007-2008 data set 
where the pinger was freely swinging. The compar- 
ison of these waveforms with those of Figure 3(b) 
shows that waveforms recorded by the same chan- 
nel are very similar even when the pinger is moved 
to a different hole. 

4.2. Pinger data analysis 

For all analyses described below, we selected data 
recorded by each channel when the pinger was at 
the same depth as the sensor, in every instrumented 
hole. This selection provides 49 combinations (four 
times three SPATS channels on each of the four 
strings, plus one HADES channel) which can be 
used for the same number of independent attenua- 
tion measurements. This data set provides the most 
systematic-free attenuation analysis, because it al- 
lows us to compare recordings of the same sensor 
as the transmitter is moved among holes at vary- 
ing distances from the sensor. Moreover, the sensor- 
transmitter configurations that were selected fea- 
ture minimal variation in both polar angle and az- 
imuthal angle of the arriving signal at each sensor, 
as the pinger is moved from hole to hole. We em- 
phasize that for each single measurement only data 
recorded by one single channel have been used. This 
means that we are not sensitive to variation of sensi- 
tivity between channels and we can neglect the sen- 
sor response function, unknown in our case, assum- 
ing the following: 

- The sensor response is constant throughout the 
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Fig. 3. Waveforms recorded by the same channel when the 
pingcr was stopped while lowering and while raising, for 
two different holes at a distance of 243 m (top) and 336 m 
(bottom). Clock drift correction and averaging of the 180 
recorded pulses have been applied. In the top plot, one of 
the waveforms has been shifted in time to overlap with the 
other one. This small shift is due to a slight difference in the 
pingcr depth between the two runs (see [21]). 



pinger data taking. This is supported by the fact 
that the noise spectra measured by the sensors 
has been demonstrated to be very stable in time. 
- The sensor output is linear in amplitude with 
respect to the input amplitude. This has been 
demonstrated in the laboratory for signals which 
are within the range of amplitudes considered 
here. Waveforms which are saturated are ex- 
cluded. 

We also make two additional important assump- 
tions, which are supported by the similarity of the 



waveforms recorded by each channel for multiple 
pingcr distances and by the analysis of the received 
signal in the frequency domain (see section 4.2.2 and 
in particular Fig. 8): 

- The pinger emission is constant throughout all 
sets of measurements taken. 

- The sound transmission in the medium is not af- 
fected by dispersion, i.e., the frequency content of 
the pulse is independent of the traveled distance. 
Two types of analyses have been performed on the 

pinger data: the energy of the full waveform calcu- 
lated in the time domain and the energy calculated 
in the frequency domain. 



4.2.1. Time domain pinger analysis 

The energy analysis integrates over the complete 
waveform. The full waveform recorded when the 
pinger was stopped at the depths of the sensors 
was processed as explained in Section 3.1. The 
high quality averaged pulse obtained was used to 
calculate the energy Es+n for each channel- hole 
combination by applying Equation 4. This includes 
both a signal and a noise contribution. The noise 
energy En is calculated using the data recorded 
immediately before and after the pingcr operation, 
which was verified to be very stable over the time. 
The processing of the waveform is done in the same 
way for both kinds of data (signal plus noise or 
noise-only). The noise was subtracted to estimate 
the signal energy E = Es+n — En- Next, the ef- 
fective amplitude was calculated (Equation 5). The 
number of valid measurements was 48 of the 49 
channels available, as one channel (HADES, at level 
190 m in string D) did not have more than two data 
points which pass the cut E > 0. 

The statistical error on the effective amplitude 
was determined by error propagation from the sta- 
tistical uncertainty of the sample amplitudes as ex- 
plained in section 3.3. To determine the systematic 
uncertainty, we calculated for each of the 48 mea- 
surements the x 2 l n d (where is the number of de- 
grees of freedom of the fit) for different hypothesized 
systematic error values and we observed how the 
distribution of these values changed. We found that 
adding 15% systematic uncertainty to each fit brings 
this distribution to have mean equal to 1 . We there- 
fore used this value as an estimate of the systematic 
uncertainties among data points for each fit, which 
can be attributed to residual azimuthal and polar 
variation of the sensor sensitivity (due to a non- 
perfect alignment of holes and stopping depths) . We 
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Fig. 4. Example fit of effective amplitude vs distance for 
sensor channel BS6-0 installed at 320 m depth. Error bars 
include the statistical and systematic uncertainty, estimated 
as explained in the text. 
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module 

Fig. 5. Attenuation coefficients (with standard error) ob- 
tained by energy analysis in time domain. The three chan- 
nels of each module are indicated by the three data points at 
and to the right of the corresponding label on the horizontal 
axis. 

also found that variation of the assumed systematic 
uncertainty within a reasonable range does not sig- 
nificantly affect the final result (see [21]). 

Taking into account the statistical and systematic 
uncertainties, we fit Equation 3 to determine the 
amplitude attenuation coefficient a. A typical fit of 
the data, from sensor channel BS6-0, is shown in 
Figure 4. 

A compilation of results obtained by fitting the 
single-channel data is shown in Figure 5. The strings 
have been sorted from nearest to farthest from the 
pinger holes and, within each string, modules have 



Fig. 6. Mean attenuation coefficient as a function of the 
depth. 

been sorted by depth. For details of each single- 
channel fit see [21]. A consideration must be made 
looking at the distribution of the points: the scat- 
ter of the data, greater than the error bars, implies 
that there are additional systematic uncertainties 
(for example, local properties of the ice, or of the in- 
terface between hole ice and sensors) which we are 
not able to identify and quantify without looking at 
the spread of the data. Nevertheless, the 48 mea- 
surements allow us to constrain the attenuation co- 
efficient within a narrow range. To take into account 
the uncertainties of the individual measurements we 
assign to each entry a weight inversely proportional 
to the error of the value: 



l 




(7) 



where a a ^ is the onc-sigma uncertainty on the value 
of a in the measurement i obtained from the fit. We 
then build a histogram using all the weighted entries 
and we take the mean and the standard deviation 
(indicated by the dashed line in Figure 5) as final 
value and error: 

(a) = 3.20 ± 0.57 km- 1 . (8) 

In Figure 6 the average attenuation coefficient is 
shown for each depth. No clear depth dependence is 
observed. 

4.2.2. Frequency domain pinger analysis 

The data selected for the frequency domain anal- 
ysis are the same as in the previous analysis, i.e. 
the attenuation coefficient is obtained from the 
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horizontal pingcr-sensor configuration. However, 
the differences in the method, including averaging 
and background handling, allow for a cross-check 
of the methodology. Moreover the analysis software 
used for the two methods was completely separate, 
allowing for an independent cross-check. 

For each channel a series of pinger pulses (I = 
1, ...,N p ; N p =180 or 144) is recorded. Each pulse 
is the sum of signal s and an additive stationary 
noise component n (uncorrelated with the signal), 
therefore the amplitude x of a sample k at time tk 
can be written as: 



ar'(tfc) = s l (t k ) + n l (t k ). 



(9) 



The discrete Fourier transform yields the complex 
numbers 
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Fig. 7. Example of frequency spectra of signal plus noise, 
pure noise and signal after noise subtraction 
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where (k — 1, ...,M) is the sample number within 
each pinger pulse (M = 20 or 25). 

The average noise, N, is estimated from the off- 
pulse portion of the waveforms (according to the 
methods described in [23-25]) by averaging over a 
set of pure noise samples taken from each waveform. 
To avoid any overlap with the signal, the noise in- 
tervals are taken before each recorded pulse. Exam- 
ples of a raw signal plus noise spectrum, a pure noise 
spectrum and a signal spectrum arc shown in Fig- 
ure 7. 

As can be seen from Figure 8, the spectral shape 
is approximately constant for the same sensor, but 
attenuated with increasing distance, i.e. for differ- 
ent pinger hole measurements. The peak at 10 kHz 
reflects a characteristic peak in the sensor response. 

We compute the average signal spectrum density 
S m after subtracting the noise density from each 
pinger pulse: 



N ^ 



N„ 



(12) 



where the sum extends over the pulses and applies in 
the region 5-30 kHz. Next we compute the effective 
amplitude A e f f via Equation 6 and apply the same 
fit procedure as in the time domain analysis. 

The weighted mean of a for the available channels 
from the frequency domain analysis is 



I 1 1 1 1 I 1 1 1 1 I 
, SPATS AS5(1) 





3.75 ± 0.61 km" 



(13) 



25 30 

frequency [kHz] 

Fig. 8. Example of typical noise-subtracted spectra. The two 
spectra are taken from the same sensor, but for different dis- 
tances. The Fourier magnitudes arc multiplied by distance, 
which then allows a direct estimate of the attenuation, as 
described by Equation 3. 

where the uncertainty is given by the width of the 
weighted distribution. The difference of the value 
from that of the time domain analysis is explained 
by the different method used to separate signal from 
noise; the two results are nevertheless compatible. 

A similar study has been done using as input the 
averaged waveforms, separating the energy contri- 
bution for the frequency intervals between 5 and 
17 kHz, between 17 and 35 kHz, and for the whole 
interval 5 to 35 kHz. We investigated the possible 
presence of a trend in the variation of a over each 
interval compared to the one obtained using the full 
spectrum, and we have not found any significant fre- 
quency dependence. In some cases the attenuation 



11 



below 20 kHz is stronger than around 30 kHz (as 
seen for example in Figure 8); in other cases the vari- 
ation is opposite. In all cases, the obtained a values 
agree within uncertainties [21]. 

Both SPATS and HADES sensors obtain the same 
shape for the frequency distribution, despite their 
different frequency responses. The recorded spec- 
trum in both cases is dominated by the pinger fre- 
quency content rather than by the response of the 
sensors. 

5. Inter-string measurement 

In addition to the measurement performed with 
the retrievable pinger operated in water, we mea- 
sured the acoustic attenuation using the frozen-in 
SPATS transmitters. Two independent analyses of a 
single data set were performed. The processing tech- 
niques used to average the waveforms and estimate 
the noise-subtracted signal energies use the same al- 
gorithms, in the time domain (see Section 3.1) but 
were implemented independently to cross-check one 
another. The signal energies determined by the two 
analyses are consistent. 

5.1. Inter-string data acquisition and processing 

Data were taken for each transmitter recorded by 
each sensor in the SPATS array. In each run, one 
transmitter was pulsed at a 25 Hz repetition rate for 
40 s. Data from all inter-string transmitter-sensor 
combinations were collected over a two-day period 
(April 1-2, 2009). The transmitters require several 
seconds to reach steady pulse-to-pulse performance, 
likely due to self-heating of the electronics during 
initial operation followed by temperature equilibra- 
tion. The initial acoustic pulse amplitude is some- 
what larger than in the steady state, followed by de- 
cay to the steady-state amplitude with a time con- 
stant of ^2 s. To be sure we are recording in the 
steady state, we start the sensor recording lis after 
the transmitter begins pulsing. In each run, all three 
channels of a single sensor module are recorded con- 
tinuously for 20 s, enough time to record 500 trans- 
mitter pulses. The three channels are sampled syn- 
chronously, at 200 kilosamplcs per second on each 
channel. 

The ambient acoustic noise conditions at the 
South Pole have been determined by SPATS to be 
very stable [26] . In particular we checked that they 
were stable during the two-day inter-string data- 



taking period. Several raw noise runs were taken on 
each channel, interspersed among the transmitter 
recordings, using the same sampling frequency and 
time duration but with no transmitter pulsing. The 
noise runs were processed with the same waveform 
averaging algorithm as the transmitter runs. From 
this data the DC offset (fi) and the standard de- 
viation of the noise samples (c) were calculated 
for each sensor channel. Each of \x and <r, on each 
sensor channel, was stable at the few-percent level 
during the two-day period. 

For each sensor recording, the mean waveform 
was processed to determine the noise-subtracted sig- 
nal energy following the algorithm described in Sec- 
tion 3.1. 

5.2. Inter-string single-depth direct analysis 

This analysis uses a single transmitter recorded by 
all sensor channels at the same depth as the trans- 
mitter. A single transmitter is used because the dif- 
ferent transmitters are known to have different in- 
herent transmittivity and perhaps different coupling 
to the ice with respect to one another. Only those 
sensor channels at the same depth as the transmit- 
ter were used in order to mitigate effects due to the 
unknown change of transmittivity with varying po- 
lar angle. 

For each transmitter there are typically three sen- 
sor modules, each with three channels, located at 
the same depth as the transmitter but on other 
strings. The acoustic attenuation was measured for 
each transmitter in the array, after applying the fol- 
lowing two quality cuts to the data. First, the sta- 
tistical uncertainty of the effective amplitude at a 
given channel was required to be 20% or smaller, 
in order to be considered a "good" channel and in- 
cluded in the fit. Second, there had to be at least one 
good channel at each of at least two distances from 
the transmitter, at the same depth as the transmit- 
ter, in order for a fit to be performed for a given 
transmitter. Note that only five of the nine instru- 
mented depths arc instrumented on all four strings. 
Twelve of the 28 SPATS transmitters met these two 
selection criteria. 

For each good sensor channel, the quantity y (de- 
fined by Equation 2) was determined following the 
procedure described in Section 4.2. We performed 
a linear fit for each of the 12 transmitters to deter- 
mine the acoustic attenuation coefficient a and the 
free normalization parameter b related to the emis- 
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Fig. 9. Compilation of attenuation fits for the 12 transmitters 
used in the single-depth inter-string analysis. The weighted 
mean is shown as a solid horizontal line, and the onc-sigma 
uncertainty of the global fit is indicated with dashed lines. 



sion strength of the transmitter. 

Both statistical and systematic uncertainties were 
calculated for each y value for each channel recording 
each transmitter. The statistical contribution was 
determined by propagating errors from the statisti- 
cal uncertainty of the effective amplitude. The sys- 
tematic uncertainty is dominated by the unknown 
relative sensitivity of the sensor channels. We esti- 
mate the channel-to-channel variation in sensitivity 
by treating the pinger results as an in situ calibra- 
tion of the sensor channels. Since the pinger analysis 
fits the b parameter for each channel independently, 
e b can be taken as a measure of the sensitivity of 
each channel. The results of this in situ calibration 
were used to estimate the systematic uncertainty 
in y for the inter-string attenuation analysis, which 
combines different sensor channels for each fit. The 
absolute systematic uncertainty in y was 0.97. This 
was added in quadrature to the statistical uncer- 
tainty of y for each channel. 

Results of the fits for the 12 transmitters are com- 
piled in Figure 9. For each transmitter i, the best fit 
attenuation coefficient a* and corresponding uncer- 
tainty Oi are shown. Many of the single-transmitter 
fits individually give a result consistent with zero at- 
tenuation. However, the 12 results can be combined 
to improve the precision of the measurement. De- 
termining the weighted mean of the measurements 
(weighting by 1 /of) gives a global best fit of 



5.3. Inter-string multi-depth ratio analysis 

None of the permanently deployed transmitters 
or sensors has been calibrated in ice. Therefore, 
both an unknown inherent sensitivity Sj and trans- 
mittivity % enter into the equation for a single 
inter-string amplitude measurement. With these 
variables, Equation 1 can be written for the combi- 
nation of transmitter i and sensor j as: 



(15) 



where Aij is the recorded amplitude of the pulse 
transmitted by transmitter i as detected by sen- 
sor j and dij is the distance between transmitter i 
and sensor j. If we then take two transmitter-sensor 
pairs, transmitters Tj and Tfc heard both by sensors 
Sj and Si, it is possible to construct a ratio [20] of 
amplitudes: 



. AijAki dijdki 
\n{n A TZ d ) = In \ -f— 



-aD x + 6, (16) 



3.16 ±1.05 km" 



(14) 



where TZa and IZd are ratios of amplitudes and dis- 
tances respectively, D x = ([dij - da] - [dkj - dki]) 
and b is a free fit parameter introduced to allow for a 
systematic shift in y. All three channels of both sen- 
sor modules are required to have A > 1.5tr s t a t(^4), 
where cr s t a t (A) is the statistical error on the ampli- 
tude A. Dead and saturated channels are excluded. 

One single measurement would yield the attenua- 
tion coefficient a if the transmitters and sensors were 
all perfectly isotropic. However, the azimuthal ori- 
entation of each sensor module is unknown and the 
transmitter signal is known to vary significantly with 
polar angle. Both effects need to be minimized and 
need to be accounted for. Therefore, for the anal- 
ysis presented here, only transmitter-sensor combi- 
nations from the lower neighboring levels were used: 
(190, 250), (250, 320), and (320, 400) m depth. This 
limits the difference in polar angle for an amplitude 
ratio to a maximum of 32 . The maximum differ- 
ence in azimuth angle is ^120°. An average sensor- 
module response is used so that the azimuth effect 
is less significant. A total uncertainty on the sin- 
gle amplitude ratio of 100% has been estimated by 
studying the spread of the b values obtained from the 
pinger analysis, as described in Section 5.2. On top 
of that, the variation due to the angular variation of 
the SPATS transmitters needs to be included. This 
was studied in the laboratory, where a maximum 
amplitude variation of 40% in transmitter emission 
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was observed for the inter-string geometry discussed 
here. For the ratios with amplitudes at minimal po- 
lar angles, the assumed uncertainty on a single ratio 
is conservative. 




D x [m] 



Fig. 10. The 2- level inter-string ratios as a function of D x . 
The spread on the points is an indication of the residual 
polar and azimuthal systematic dependencies. The solid line 
shows the best fit; the dotted line shows y = for reference. 

Figure 10 shows all 172 ratios with the estimated 
systematic error bars and the best fit obtained by 
constraining the offset (i.e. the parameter b of Equa- 
tion 16) to be 0. As mentioned above, the figure re- 
flects the fact that the systematic errors are overes- 
timated. However this does not influence the central 
value of the linear fit coefficient a. The fit results in 
a value 

a = 4.77 ± 0.67 km- 1 . (17) 
6. Transients measurement 

In addition to the signals from the retrievable 
pinger and the frozen-in transmitters, we can use 
ambient transient events to estimate the attenuation 
coefficient. SPATS runs in a transient data acquisi- 
tion mode during most of each hour. Three channels 
from each of the four strings are monitored for large- 
amplitude events. A simple 5cr threshold trigger is 
used, where a is the Gaussian noise amplitude. The 
Gaussian noise is determined on each channel and 
accordingly a different absolute threshold is used on 
each channel. Events above thresholds are recorded 



with an absolute GPS time stamp; this occurs with 
a frequency of ~1 Hz per channel. The hits from the 
four strings (where a hit is defined to be a times- 
tamped waveform recorded from an individual chan- 
nel) are time ordered offline and processed through 
an algorithm to find hits within a coincidence time 
window of 200 ms (the time necessary for a pres- 
sure wave signal to cross the SPATS array). Each 
cluster of hits occurring within the coincidence time 
window is an event. For events with more than four 
hits and at least three strings, a vertex location is 
calculated. 







- , , , 1 t t t I t 1 1 


1 t t t I i i i 


1 2 

sensor channel C 


3 4 J 
Sb-1 


~ , , , I i i i I i i i I i i i I i i i 


sensor channel B 


S5-1 




- , , , l. i i i l i i i 





2 3 
sensor charnel AS5-1 



Fig. 11. Waveforms produced by sound from rcfrcczing ice in 
Hole 19 received at the 250 m sensors of Strings C, B and A 
at a distance of 336 m, 632 m, and 744 m respectively. The 
time scale refers to the trigger time which is independent for 
each channel. 

Using this data taking mode for 45 minutes out of 
every hour, data have been collected since the end 
of August 2008. Acoustic signals have been observed 
for several days from re- freezing IceCube holes with 
precisely known x-y coordinates. Because some of 
these holes were also used for collecting pinger data, 
the corresponding sensor calibration constants cal- 
culated from those data can be applied to signals 
whose source is localized to near the correspond- 
ing pinger stop depths. Thirteen appropriate events 
have been found with source location at the co- 
ordinates of holes 19, 20 and 28, and depths be- 
tween 230 and 270 m, which fit in the dynamic range 
of the SPATS sensors at all strings involved. The 
waveforms of an example event are shown in Fig- 
ure 11. Single-event analysis using effective ampli- 
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tudcs derived from energy in the time domain (see 
Section 4.2) leads to a mean attenuation coefficient 

a- = 3.64 ± 0.29 km" 1 , (18) 

compatible with the results from both the retriev- 
able pinger and the frozen-in transmitters. The rela- 
tively small errors in comparison to the pinger anal- 
yses may be due to the fact that here only sensors 
installed at one depth (250 m) have been used. The 
error shown is the statistical error only and neglects 
any additional uncertainty of the sensor calibration 
in comparison to the pinger analysis. The frequency 
spectrum of the transient events needs further study 
but contains contributions up to 80 kHz. 

7. Discussion and outlook 

The results from the various analyses described 
above are summarized in Table 2. 

Tabic 2 

Summary of results from various attenuation analyses. 



Analysis Measurements a [km 



Pinger (time dom.) 


48 


3.20 


± 


0.57 


Pinger (freq. dom.) 


39 


3.75 


± 


0.61 


Inter-string direct 


12 


3.16 


± 


1.05 


Inter-string ratio 


1 


4.77 


± 


0.67 


Transients 


13 


3.64 


± 


0.29 



All of the analyses yield consistent results. How- 
ever, we cannot simply average them, due to the 
quite different systematic effects (as well as corre- 
lations among them) which affect each of them as 
discussed in detail in previous sections. In Table 3 
the different effects are summarized and attributed 
to the individual studies they affect. The pinger ap- 
proach uses the same sensor and the same transmit- 
ter at the same zenith angle and nearly the same 
azimuthal angle at different distances, thereby mit- 
igating all effects connected with unknown sensi- 
tivities of emitter and receiver. For this reason we 
believe it provides the most reliable result. Relying 
on the same data set, the pinger analyses in time 
and frequency domain are not independent but cross 
check one another using different signal and back- 
ground averaging and subtraction methods. Because 
there is a higher number of measurements passing 
selection criteria in the time domain analysis than 
in the frequency domain analysis (see Table 2), we 



quote this result as our best estimate of the attenu- 
ation coefficient: 

(a) = 3.20 ±0.57 km" 1 . (19) 

which expressed as an attenuation length is: 

(A) = 312tg m. (20) 

The results of all analyses are consistent with an 
attenuation length of ~300 m ± 20%. 

Up to 30 kHz no strong frequency dependence of 
a has been found. There are also no indications of 
depth dependence of the attenuation up to 500 m 
depth. 

Our measured value for the attenuation coeffi- 
cient is an order of magnitude larger than expected. 
In [14], it was estimated that South Pole ice grains 
are sufficiently small that Rayleigh scattering is neg- 
ligible and attenuation is dominated by absorption 
due to proton reorientation. However, new data [27] 
from the SPRESSO site near the South Pole indicate 
that the ice grains are larger than previously esti- 
mated. Because Rayleigh scattering increases with 
the cubic power of the grain length, it could be that 
the attenuation we have measured is dominated by 
scattering. Another possible mechanism of absorp- 
tion, not previously considered, is internal friction 
at linear crystallographic dislocations. 

The weak frequency dependence observed in our 
measurements below 30 kHz disfavors the Rayleigh 
scattering hypothesis. New pinger measurements 
have been taken in the 2009-2010 South Pole sum- 
mer in order to clarify the situation. A modified 
pinger, emitting a sequence of acoustic pulses at 
three different frequencies up to 60 kHz, was de- 
ployed up to 1000 m depth in three additional bore 
holes. Analysis of these data is currently under- 
way, and should allow for a more conclusive study 
of both the frequency and depth dependence of 
acoustic attenuation in South Pole ice. 

Given some of the inherent advantages of the 
acoustic technique relative to the radio technique 
(such better shielding from anthropogenic surface 
backgrounds), this study was undertaken primarily 
to determine whether the acoustic method could be 
a basic ingredient of a 100 km 3 hybrid detector for 
ultra-high energy neutrino detection at the South 
Pole. The design chosen in [15] with horizontal 
string distances of order 1 km will not reach the 
necessary sensitivity, given the 300 m attenuation 
length reported here. New geometries are under 
study to assess whether detectors with closer string 
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Table 3 

Systematic effects present in individual analyses. 
Systematic Effect 

Channel-to-channel sensitivity variation no yes no 

Azimuthal sensitivity variation minimal yes yes 

Polar sensitivity variation minimal no yes 

Channel-to-channel transmittivity variation 
Azimuthal transmittivity variation 
Polar transmittivity variation 



spacing but larger area at smaller depth could be 
an alternative solution. Measurements are also in 
preparation to measure the absolute noise level of 
the Gaussian noise, which is the most important 
remaining ingredient to determine the feasibility 
of the acoustic method for detecting neutrinos in 
South Pole ice. 
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pinger inter-string inter-string transients 

single level ratio 

minimal 
minimal 
minimal 



no no no no 

no yes yes no 

no no yes no 
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x10 3 
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